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In this communication, we would like to comment on 
sonic results in a recent Letter by Kafesaki et al. 

In the authors stated "Up to now the extensive 
studies of bubbly water (or other bubbly liquids) have 
been analyzed in the framework of single scattering" . 
In fact, acoustic scattering by bubbles has also been 
studied in terms of multiple scattering both extensively 
and intensively in the literature. This includes both 
theoretical (e. g. [p|-p^) and experimental works (e. g. 
|l|,|l|). The work of Fold y f^ , Lax Waterman 
et al. [0, and Twersky [[18[ serves as a solid foun- 
dation to the multiple scattering of acoustic waves by 
many scatterers in general. The textbook by Ishimaru 
jigf l presents an excellent account of various theories 
about the multiple scattering. The authors further sug- 
gested "Thus, acoustic waves in bubbly liquids (espe- 
cially water), in addition to being very important per 
se, provide an almost ideal case for examining in detail 
the localization question". We note that the research 
on multiple-scattering induced acoustic localization in 
random bubbly water is rich (e. g. |^,|,0|o|-|2|] ) . It 
has been shown that acoustic waves can be 

localized within a range of frequencies slightly above 
the natural frequency of an air-bubble in water, and 
the localization behavior is insensitive to the configu- 
ration of bubble clouds . When waves are localized 
in bubbly water, a global coherent behavior appears. 
This is a distinct feature differentiating the localization 
effect from residual absorption |^,^,^ . The localiza- 
tion disappears when the air-bubble volume fraction is 
lower than about 10^^, in agreement with an earlier 
prediction Q]. 

To the end of their Letter, Kafesaki ||l[ stated "The 
most significant and novel results we obtained are 
shown in the panel of Fig. 3(b)...". The results shown 
in their Fig. 3(b) are a natural extension of the re- 
sults of |11|. The formulation in [ll| stems from the 
work of Foldy jl^] and Twersky and has been doc- 
umented in detail in | |25| , |28[ |. The approach has also 
been compared favorably with others in certain situa- 
tions [P5|,P8|-30|, and the results agree with the physical 



intuition (e. g. [pl|). Though not exact, the approach is 
applicable for most frequencies considered in ^ . There 
will be a small inaccuracy when urs/c exceeds 0.45 and 
significant inaccuracy is expected when Lors/c > 1. 

Using the formulations in fl^,^ for acoustic trans- 
mission and using the formulation in |32| for comput- 
ing the acoustic band structure, we examine the results 
in Fig. 3 of |l[]. It should be noted that the formula- 
tion in Iprll only takes into account the pulsating mode, 
while the formulation in incorporates all possible 



vibrational modes of the bubbles. The two approaches 
agree well except at higher resonance frequencies of 
ujvs/c « 0.48 and 0.77 [§. Thus, in the following, the 
approach in will be used. 

Employing the same parameters as in [Q, we re- 
produce partial results in Fig. 1. For comparison, we 
show the result for a single random configuration (a) 
and as well as the ensemble averaged result (b). It 
is clear that whether making ensemble average makes 
not much difference for Lorg/co < 0.6. In particular, 
the ensemble average basically plays no role for fre- 
quencies located around the complete band gap. Com- 
pared to Fig. 3 in [0, there are several differences. (1) 
The complete band gap is a little wider in the present 
case, but is consistent with a transmission calculation, 
following |22j. To further inspect the difference, we 
have also computed the band structures at the air void 
fraction of 10%, yielding nearly the same answer as 
in Fig. 1(a) of ||l|. (2) In contrast to |jl|, the trans- 
mission in the present computation changes sharply 
at the frequencies slightly above the natural frequency 
of an air- bubble in water [Q, which in fact leads to 
a first order phase transition from the extended state 
to the localized state Exploring the difference, 

we have considered possible effects of the locations 
of transmitting source and receiver on the transmis- 
sion, and find that the transmission is insensitive to 
the positions of the transmitter and the receiver for 
frequencies between ujrs/co — 0.0136 and about 0.3, 
indicating that the discrepancy in the transmission be- 
havior between the present results and that from ^ 
is not caused by a difference in the source or receiver 
location. The transmission at frequencies higher than 
roughly wr^/co = 0.3, however, is indeed sensitive to 
the location of either transmitter or receiver. (3) In 
Fig. 3 of §, the range of inhibited transmission coin- 
cides with the band gap, while in our case this is not 
the case. In addition, our results show that the ran- 
domness tends to smear out the transmission peak near 
the upper band edge, consistent with 2D cases . At 
first sight, the present inconsistency between the range 
of the band gap and the range of the inhibited trans- 
mission seems to suggest errors in our computation. 
But, in fact, the the disagreement is due to the effect 
of the finite sample size. To elaborate. Fig. 2 shows the 
results for different sizes of bubble clouds, where the 
effect of finite sample size is obvious. The range of the 
transmission inhibition shall agree with the band gap 
when the sample size tends infinity. 

When the multiple scattering is switched off, the 
effect of ensemble average becomes more prominent. 



comparing the dotted lines in the right panels in 
Fig. 1(a) and (b). Again we see discrepancies with 
Fig. 3 in For example, the present result for the 
lattice case starts to flat out at about wr^/co = 0.3, 
while that in [Q at about 0.4. The randomness tends 
to broad the transmission peak located around the 
natural frequency of the breathing mode of the air- 
bubble. The two transmission peaks at wrs/co — 0.48 
and 0.77 shown in are absent from the right panel of 
Fig. 1. This is because the results in Fig. 1 did not take 
into consideration higher vibrational modes. When the 
higher modes are considered using |28| , the overall re- 
sults in the right panel of Fig. 1 will not change much, 
except that two sharp transmission peaks will appear 
as in Fig. 3(c) of §. 

Discussion with Dr. Pigang Luan is greatly appre- 
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Figure 1 (a) Left panel. Acoustic bands in a sim- 
ple cubic array of air-bubbles in water, rg is the 
bubble radius, u is the angular frequency and cq 
is the water velocity in the water. Air volume 
fraction is 1%. Middle panel. Reduced trans- 
mission coefficient through the bubble array in 
ordered (solid line) or one random (dotted line) 
placement. Right panel. The same as in the 
middle panel, but with the multiple scattering 
switched off. (b) The same as in (a), but the 
transmission in the case of random placement is 
averaged over 100 random configurations. 

Figure 2 The same as in Fig. 1, but for various sizes 
of bubble clouds in ordered placement (a) and as 
well as for one random arrangement (b). 
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